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We analyze a model that interpolates between scale-free and Erdos-Renyi networks. The model 
introduced generates a one-parameter family of networks and allows to analyze the role of structural 
heterogeneity. Analytical calculations are compared with extensive numerical simulations in order 
to describe the transition between these two important classes of networks. Finally, an application 
of the proposed model to the study of the percolation transition is presented. 



PACS numbers: 



i.75.Fb, 05.70.Jk 



I. INTRODUCTION 

In the last several years, graph theory has experienced 
a burst of activity as many real systems can be repre- 
sented and modeled as networks [l|, Q, Q . A network is 
made up of vertices representing the interacting elements 
of the system and of edges that stand for the interactions 
among them. Network modeling comprises the analysis 
and characterization of the structure of networks as well 
as their modeling in terms of generic models aimed at 
reproducing the features found in real systems |1, 13, 13 ■ 
The second important branch has to do with dynamics 
on networks. This has lately attracted the attention of 
many scientists as it is ultimately related with the func- 
tioning of the system that is being modeled |3, 3 ■ 

In this paper, we deal with the first of these areas 
of research. The seminal paper by Barabasi and Albert 
[3, 1^ , showed that many real world networks can not be 
described by graphs where the connectivity distribution 
(i.e., the probability that a given node has a given num- 
ber of links) follows a Poisson-like distribution. Indeed, 
Barabasi and Albert showed that most real networks are 
heterogeneous in the sense that the probability that a 
node is connected to k other nodes follows a power-law 
distribution P{k) ~ k~'^, where 7 usually lies between 2 
and 3. These networks were termed scale-free networks 

QiiiQ. 

Soon afterwards, many studies have dealt with the 
analysis and characterization of models that generate 
scale-free networks, along with other global and local 
topological properties found in real networked systems 
[3,|3,HSIl3- ^^ particular, models based on the mech- 
anism of preferential attachment (PA), no matter if the 
network is growing or not, have been extensively studied 
in the last years. There are some models in which the 
PA rule is limited to a neighborhood due to g eographic 
constraints [ll| or lack of global knowledge [l3| , or where 
its linear character is investigated [l^- While today we 
have recognized that preferential attachment is not a nec- 



essary condition for the formation of scale-free networks 
[lj|) it seems to be clear that it is an important mech- 
anism. Indeed, most of the existing models intrinsically 
incorporate a preferential attachment like rule. On the 
other hand, uniform random linking of nodes on growing 
networks gives rise to networks where the degree distribu- 
tion decays exponentially fast with the degree k, thus pro- 
ducing homogeneous networks with a well defined (and 
meaningful) average value for k \iaL flq . 

The combination of the two rules, i.e, uniform and pref- 
erential linking, have been also analyzed in several mod- 
els for interpolating between scale-free and exponential 
networks rq. For instance, Liu et al 1»| have studied a 
model in which the probability of establishing new links 
goes as a linear combination of both in such a way that 
a new link is established between a node i and a new 
one proportionally to {1 — p)ki + p, where p weights the 
contribution of the two mechanisms. However, in previ- 
ous models of this sort, there is an assumption that does 
not apply always. It has to do with the fact that the 
network always grows around a single component of con- 
nected nodes and uniform or preferential links from the 
emerging nodes are always made with elements belong- 
ing to this unique cluster. This single component grows 
linearly in time until it reaches the size of the network. 
Since there are no clusters of nodes other than the giant 
component, the models can not account for phenomena 
such as the coalescence of small networks into a larger 
one, nor for situations in which more than one node is 
added to a preexisting structure at each time step, fea- 
tures that may be relevant in social, economic and other 
networked systems. 

In this paper, we analyze a model that interpolates 
between Erdos-Renyi (ER) and scale-free (SF) networks 
as far as the degree distribution is concerned through a 
tunable parameter. By construction, new links are not 
always established with nodes prcviouly incorporated to 
the network. We explore analytically and numerically 
the time behavior of nodes attachment as well as of the 



degree evolution. We find that, depending on the inter- 
play between uniform and preferential linking, the tran- 
sition from an ER like network to an SF one is smooth or 
more abrupt. We finally discuss other topological prop- 
erties and perform a numerical percolation study that 
highlights the differences in their structure. The model 
presented here is useful as it provides a unique recipe 
to go progressively from homogeneous to heterogeneous 
topologies as well as for exploring the interplay between 
them. 



II. THE MODEL 

The model introduced in this work generates a one- 
parameter family of complex networks. This parameter, 
a e [0,1], measures the degree of heterogeneity of the 
final networks. Let us assume the final size of the network 
to be ri. The network is generated in the following way: 

(i) Start from a fully connected network of jtiq nodes 
and a set U{Q) of (fi — mo) unconnected nodes. 

(ii) At each time step choose a new node j from U{Q). 

(Hi) This node makes a link in two ways: 

(a) With probability a it links to any other node 
i of the whole set of O — 1 nodes with uniform 
probability 



a ^ 



n: 



uniform 



= {n-i)-' 



(1) 



(b) With probability (1— a) establish a link follow- 
ing a preferential attachment strategy, that is, 
the probability for any other node i to attach 
to node j is a function of its connectivity as. 



n 



PA 



F(A:,; 



(2) 



where different choices for the functional form 
of F(a;) are analyzed below. 

(iv) Repeat m times step (Hi) for the same node j . 

(v) Repeat U{Q) = {Vt — mo) times steps (ii) to (iv). 

A schematic plot of the linking procedure at step (Hi) is 
shown in Fig. ^ The above rules allow for the coexis- 
tence of two classes of nodes. On one hand, there are 
nodes with at least one link. This set will be referred to 
henceforth as the connected set N{t) J^]. On the other 
hand, there is another set U{t) of isolated nodes such 
that its size is i7 — N{t). At variance with other models 
in which there are only nodes with connectivity different 
from zero and thus the connected component grows lin- 
early with time, the above rules allows the addition of 
more than one node to the set N{t) as a result of ran- 
dom linking. Therefore, we expect the time dependency 
of N{t) to be highly non trivial. 
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FIG. 1: (Color online). Schematic representation of the gen- 
eral procedure for generating the networks. With probability 
a one of the m links can be made with any of the nodes (and 
with the same uniform probability) that will take part in the 
final network. On the other hand, with probability (1 — a) 
the link will be made only with those nodes that form the 
connected set at that time because the node will choose a 
preferential linking strategy. 



III. NETWORK GROWTH AND DEGREE 
EVOLUTION 



In order to describe the evolution of the nodes degree, 
one has to consider the functional form of F(x) for the 
preferential attachment probability Q. However, we can 
take into account some previous considerations that do 
not depend on the particular form of F(x). 

First of all, it is useful to consider two kind of links in 
order to analyze the model. Namely, the ones that arise 
from a uniform random choice, fc", and the remaining, 
fcP'^, corresponding to the implementation of the prefer- 
ential attachment rule. The dynamics of fc" is completely 
independent of the dynamics of the PA links, kP^, but 
the opposite is not necessarily true. From this, it follows 
that the probability that one node has k^ uniform links, 
P^{k^), is a Poisson distribution with (fc") = 2am. 



F"(fc") 



{2am)'' 



A:"! 



(3) 



As a consequence, we will concentrate on analyzing the 
growth dynamics of the PA links for the studied models. 
It is particularly interesting to study at this point how 
uniform random linking affects the evolution of the con- 
nected set since this is completely independent on the 
specific PA rule considered. This feature represents one 
of the main differences between the studied model and 
other previous mechanisms used to generate growing net- 
works [3,yjl3- That is, in our model nodes are not incor- 
porated to the connected set at a constant rate (like e.g. 
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FIG. 2: (Color online). Size of the connected set N{t) as 
a function of time. Solid lines correspond to the analytical 
results (eq. ^) and points are the Monte Carlo results of 
network construction (employing model A (sec. 1111 AH . The 
comparison is made for Q = 10^ and several values of a. The 
parameters of the model are set to A — m — mo — 1. 



in the standard Barabasi- Albert model) due to the possi- 
bility of adding new nodes from U{t) by applying uniform 
linking at time t and therefore the set U{t) ^ ^(0) — t. 
We can easily derive the evolution of the connected set 
size, N{t) = n — U{t), for any value of the parameter a. 
For this, we consider the growth of the connected set at 
each time step, i.e. when a new node oiU(0) throws its 
m links 



N{t+l) = N{t) + 



n-N{t) 
n- (t + mo) 



0-^)^<^' 



In the above equation the second term on the right ac- 
counts for the probability that the new node (which is 
throwing its m links) of U{Q) does not belong already 
to the connected set at time t (due to the possible uni- 
form links obtained from previous nodes of iY(0) already 
connected to the connected set N{t)). Besides, the third 
term on the right describes the probability that any uni- 
form link thrown by the node is directed to a node be- 
longing to U{t). These two terms account for the growth 
rate of the connected set. We can consider that both time 
and N(t) are continuous variables and make the time step 
small enough in order to obtain the corresponding ODE 
associated to cq. Q, whose solution is given by 



N{t) ^n + it + mo~ n)e-""''/^ 



(5) 



The agreement between this calculation and Monte Carlo 
simulations is shown in Fig. |21for different values of a and 
a preferential attachment as described in what follows 
(model A) IIII Al It is worth noting the highly nonlinear 
behavior of N(t), at variance with models in which its 
size changes at a constant rate. 
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FIG. 3: (Color online). Model A. Monte Carlo simulation 
(points) versus mean field (lines) results for k^^{t — fl) as a, 
function of the birth time to for different values of a. The 
parameters of the model were Q = 10^ and A = m = mo = 1. 
The statistics of the Monte Carlo simulations were performed 
using 10* networks for each value of a. 



We formulate below two different ways to implement 
the preferential attachment rule, which give rise to dif- 
ferent behaviors. In both models we will consider that 
the PA probability of a node j depends only on the PA 
links of the node, k^^. This new separation between PA 
links and uniform ones introduces a higher differentiation 
between the two simultaneous kinds of link dynamics im- 
plemented here allowing us to manipulate (as shown be- 
low) the degree of correlation between them. The two 
models interpolate between scale-free and Erdos-Renyi 
topologies but the structural transition is quite different 
(as we will show in section IVIl. 



A. MODEL A 

In this first model we shall study a preferential at- 
tachment rule strongly correlated with the simultaneous 
uniform random linking. First, we consider that the PA 
probability of a node i is linear with the incoming PA 
degree of the node, kf^, that is, those links received by i 
when other node launches (in average) {l — a)rn links fol- 
lowing the PA rule. This particularity of the PA rule was 
already considered by Dorogovtsev et al ^^. Besides, we 
consider that when a node is introduced in the connected 
component (because cither it is chosen at random by any 
node or it is launching its m outgoing links over the rest 
of nodes) it has an initial attractiveness (or fitness) A. 
In other words, each node has an associated parameter 
Ai that is zero if the node i is not in the connected set 
and is Ai = Aifi is linked to other nodes {i.e., it belongs 



to N{t)). We further consider that the attractiveness Ai 
enters hnearly in the preferential hnking probabihty of 
node i. With these two ingredients, the expression for 
n^^ is given by 



n 



PA 



kr+A, 



(6) 



The introduction of the fitness A correlates the PA rule 
with the uniform linking in the sense that the more links 
are established uniformly (the higher a), the more 
new nodes with kf^ = are incorporated to the con- 
nected set from U{t) and hence (by the presence of A in 
the PA probability) the more candidates to obtain PA 
links are available. This can be observed from the evolu- 
tion of the connected set N{t), when a is high there are 
a lot of nodes added into N{t) at the early stage of the 
network construction so that the potential growth of the 
PA degree of the former members of the connected set 
is strongly weakened. In order to confirm these heuristic 
considerations we derive the mean field evolution for the 
incoming PA degree of a node i, kf^ 



dkf 



di 



(1 — a)m 



fcf 



A 



(1 -a)mt + AN{t) ' 



(7) 



(with the initial condition fcf^(io) = ^)- Obviously, in the 
limit a = we recover the mean field equation for the 
Generalized Dorogovtsev model J20| (which, when A = 
m, describes the Barabasi- Albert model). For a 7^ the 
influence of the uniform random linking is evident from 
the presence of N{t). The number of nodes that start 
to have the above dynamics at some time to is dN{t)/dt 
evaluated at time t — to which for a ^ is not constant 
as we have seen in the previous calculation of N{t). The 
solution of O is then given by 



kr{t = n) 



-1+exp 



(1 — a)m 



dt 



B. MODEL B 

In the second proposal the two different linking pro- 
cesses are completely independent. For this, we consider 
that nf ^ is a linear function of the (incoming and outgo- 
ing) links that appear as a product of the application of 
the PA rule. Then, kf^ will be zero until it launches its 
am PA links over the rest of the nodes, i.e. regardless 
of kf. Then, the mean field equation for the evolution of 
kf^ is given by 



(1 -a)mt + AN(t) 

(8) 
in order to obtain 

= f2) -|- am) as a 



We have solved numerically eq. © 

krit = n) {or krit = n) = knt ^ 

function of ^q. This function, along with the number of 
nodes which are incorporated to the connected set at time 
ig = to, gives the degree distribution of the PA links. We 
have compared the results given by eq. (|HJl for different 
values of a with the corresponding ones obtained by per- 
forming Monte Carlo simulations of the model (averaging 
over 10** networks for each value of a). The results, plot- 
ted in Fig. 13 show a very good agreement for the mean 
field model and the numerical network construction. As 
expected, the sooner a node is incorporated to the con- 
nected set the higher its final PA degree. However, as 
discussed above, one can observe that this gain of the 
oldest nodes becomes less important when the value of 
a grows due to the combination of two effects: (i) the 
application of the PA rule becomes less frequent and (ii) 
the fast growth of the connected set tends to make more 
homogeneous the PA probability of the nodes. 



dfcf 



dt 



(1 — a)m. 



2(1 — a)mt + Too 



(9) 



with the initial condition fc|'^(to) = {l — a)m and tg being 
the time when node i launches its to links. Solving the 
above equation yields 



fcP^(i) === (1 _ a)m 



1/2 



(10) 



Because the nodes launch their links at a constant rate 
(one node per time step), it is easy to obtain the degree 
distribution P{k^'^) 



P(fcP^) = 2(l-a)2m2(fcP^)- 



(11) 



which is simply a power law distribution with a Barabasi- 
Albert exponent regardless of the value of a. On the 
other hand, the relative weight of the power law with 
respect to the Poisson distribution in the total degree 
distribution P(fc) will be obviously affected by a (as the 
prcfactor in the above equation suggests). 



IV. NETWORK PROPERTIES 

. In this section we discuss the transition from SF to 
ER networks in terms of the global topological features 
of the networks. We have performed Monte Carlo sim- 
ulations of the two models and compared how the rele- 
vant topological measures evolve as a function of a. We 
are interested in obtaining how the different correlations 
between the uniform and PA linking rules affect several 
structural measures. To do this, we have studied the 
behavior of three magnitudes that behave very different 
in the two known limiting cases (SF and ER networks), 
namely: the degree distribution P{k), the average short- 
est path length {L) and the second moment of the degree 
distribution (fc^). 

Degree distribution - The degree distribution evolution 
is clearly different for the two models. In Fig. 0]we have 
plotted the degree distribution and the rank-degree rela- 
tion for both models. The rank-degree relation provides 
a useful tool for analysing the degree heterogeneity of the 
networks [3 and thus it is helpful when looking at the 
transition between ER and SF networks. As can be ob- 
served from Figs. EJa) and^c) the correlated model A 
shows a smooth transition from the power law (a = 0) 
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FIG. 4: (Color online). Monte Carlo results for the degree distribution P{k) and rank-degree relation for several values of a. 
(a) and (c) show the results for model A revealing a progressive increase of the tails decaying rate when q ^ 1. The results 
for model B ((b) and (d)) show how the decaying rate is not affected by a. The networks were generated with the following 
parameters fi = lO'' and m = mo = 3 (A = 3 for model A). 



to the Poisson distribution (a — 1). The decay of the 
tails (fc >> 1) of the degree distribution and the rank- 
degree relation becomes progressively faster as a grows 
revealing the decrease of the exponent of PP'^{k^'^) as 
expected from the results obtained by the analytical in- 
sights developed for model A. For model B the transition 
is completely different as it is shown in Figs. E^b) and 
0Jd). In both representations the decaying rate of the 
tails is independent of a and the transition to the Pois- 
son distribution is much more apparent for low values of 
k. In this sense one can conclude that highly connected 
nodes persist along the transition of model B while for 
model A the heterogeneity is progressively lost. 

Average shortest path length - The different evolution of 
the degree distributions observed above suggests to look 
at how the average shortest path length behaves along 
the two paths of interpolation. It is well known that 
the existence of high degree nodes makes the network 
more compact due to the possibility of finding shortcuts 
between nodes going through the hubs. Hence, the per- 
sistence of highly connected nodes determines the small 
diameter of the scale-free network. The results obtained 
are shown in Fig. EJa). As expected, the average short- 



est path length as a function of a grows slower for model 
B because the probability of finding hubs is higher than 
for the networks generated using model A for the same 
value of a. 

Second moment of P{k) - In order to obtain a quantita- 
tive measure of the evolution of the degree heterogeneity 
for the two models it is convenient to measure the second 
moment of the degree distribution, (fc^). This magnitude 
diverges (in the thermodynamic limit fl — > oo) for scale- 
free networks with exponents between 2 and 3. So, we 
expect a decrease of the heterogeneity on the path to 
ER graphs. As can be observed from Fig. Ob), model 
A shows a faster decrease of (fc^) as expected from the 
study of the degree distribution while for model B the 
transition is much smoother revealing again the persis- 
tence of highly connected nodes along the path to the ER 
limit. 

As for other properties like the clustering coefficient 
and degree-degree correlations we have checked that they 
remain unchanged irrespective of the value a and wheter 
model A or B is implemented. 
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FIG. 5: Average path length (a) and second moment of the 
degree distribution (b) as a function of a. Both quantities 
are represented normalized by their respective values in the 
ER limit. The results clearly manifest the two different tran- 
sitions of the models regarding the heterogeneity evolution 
along the interpolating path. The averaged networks had the 
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V. PERCOLATION DYNAMICS 

One of the most important differences between ER 
graphs and SF networks is given by the radically dif- 
ferent behaviors of dynamical processes that take place 
on top of them jH |22, 23, 24, 25, 26). For instance, 
epidemic spreading processes show a natural threshold 
below which the epidemic can not spread for ER graphs, 
while this threshold is absent in the thermodynamic 
limit in SF networks with a diverging second moment 
|24..25. . .26j . This kind of behavior is precisely what makes 
scale-free networks so special. 

We have implemented a percolation process on top of 
the networks generated by models A and B. It is aimed 
at simulating the random failures of a fraction / of nodes 
|2ll l22 | . By computing the size of the giant connected 
component, one can characterize the percolation transi- 
tion. Here, however, we are not interested in the transi- 
tion per se, but on the influence of the topological fea- 
tures unraveled in the previous section on the size of the 
giant component of the network. As expected, the be- 
havior of the two models is different when the limit of 
ER graphs is approached for the same values of a. 

As usual, we have analyzed the evolution of the size 
of the giant component of the network after a fraction / 
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FIG. 6: (Color online). Ratio between the size of the giant 
connected component in networks generated using models A 
(upper) and B (bottom) and ER networks when a fraction / 
of nodes is randomly deleted for different values of a. Note 
that although the transition points for models A and B are 
locate close to each other, the dependence with a is quite 
different. Network parameters are those of Fig|S] 



of the nodes (and hence their links) are removed from 
it. In this sense, it is relevant to study the relation 
G'(/)/Ger(/), where G(/) (Ger(/)) is the size of the 
giant component of the network for a j^ 1 {a — 1), 
since this will clearly unravel the two different approaches 
from the scale-free limit to the Erdos-Renyi network. As 
can be observed from Fig. |H|the differences between the 
SF and ER networks are relevant when / is high due to 
the different critical behaviors near the transition point 
[21I |23 | . More interesting for our concerns is the evolu- 
tion of the magnitude G(/)/Ger(/) as a function of a. 
For high values of / we can clearly appreaciate the fast 
approach from SF to ER supplied by model A (upper 
plot of Fig. IHJ). In fact, the size of the giant component 
is very similar to that corresponding to the ER graph for 
(a > 0.7). In contrast, the transition exhibited for the 
same values of / when model B is implemented reveals 
again a slower transition between SF and ER networks 
(Fig. El down). In particular one can observe that the 



values of G(/)/G'er(/) for a ~ 0.9 are similar to those 
for a — 0.5 in model A. 

The two different behaviors observed can be explained 
in terms of the topological measures showed in the previ- 
ous section, in particular with the behavior of the second 
moment of the degree distribution (fc^) which is known 
to play a key role for understanding the different transi- 
tions in SF and ER networks. It is known that the critical 
value of / relates to the moments of the distribution as 
l-fc = 1/ko, whereto = {{k^)/{k))-l- From Fig. Elb), 
one can see that for model A, the second moment of the 
degree distribution approaches the value obtained in ER 
limit for intermediate values of a. This is no anymore 
the case for model B, in which the ratio (fc^}/(fc^)ER ap- 
proaches 1 only for values close to a = 1. Given that the 
average connectivity (fc) is the same for all values of a in 
both models, we should expect that the behavior of a per- 
colation process near the critical point would be nearly 
the same to that of ER networks when (fc^)/(fc^)ER -^ 1- 
This is the situation for a broader range of a values in 
model A as suggested by Fig. |3b). In other words, the 
different behaviors for the two networks families observed 
in Fig. can be explained in terms of the speed at which 
(fc^) convergences to its corresponding value in the ER 
limit. 

The implemented percolation dynamics serves us to il- 
lustrate how the interpolating model presented in this 
work can be a useful tool for discovering what topologi- 
cal features are fundamental to explain the different be- 
haviors observed when more complex dynamics are stud- 
ied. Moreover, real networks are not purely scale- free nor 
completely ER, so that their behavior in front of perco- 
lation like processes lies in between these extreme cases 
as clearly shown in Fig. 



combination of uniform and preferential linking allows us 
to explore the whole path between the two limiting cases. 
An important feature of our model is that the size of the 
connected set does not grow linearly with the number of 
nodes attached to the network. This is a result of the 
novel ingredient of the model, that through the uniform 
random linking allows every node to take part on the 
network irrespective of their connectivities (even if they 
are not connected at all). We have analyzed two dif- 
ferent variants for the interpolation between the Erdos- 
Renyi and scale-free limits. On one model the transition 
is smooth, while for the other it becomes sharper. The 
analytical insights together with numerical simulations 
supported that the differences in the versions analyzed 
are rooted at the interplay between uniform and prefer- 
ential attachment. Finally, simulations of a percolation 
process have illustrated the differences in both formula- 
tions and their associated transitions (which are easily 
explained when looking at the topological properties of 
the networks). 

As for future works, the present model provides a use- 
ful tool to study the influence of the degree of hetero- 
geneity in dynamical processes of different kinds just as 
the Watts-Strogatz model have proved to do so in the 
transition from regular to random structures. In partic- 
ular, there exist open questions in phenomena such as 
the synchronization of coupled oscillators | .27j where this 
kind of model could be particularly relevant to explore 
the system's behavior in the region where homogeneous 
and heterogeneous architectures coexist. 
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